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Using the SED to locate the 7-ray emission site of powerful blazars 
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The location of the Gamma-ray emission of powerful blazars is a matter of active debate. Is the 
location within the UV emitting sub-pc scale broad line region, or farther out at pc scales where the 
molecular torus IR emission dominates? We present a diagnostic that connects three observables, 
the synchrotron and external Compton peak frequencies and the Compton dominance (the ratio 
of Compton to synchrotron luminosity) to the seed photon energy and energy density. We discuss 
encouraging preliminary results and discuss how to use our diagnostic to understand the location 
of the Gamma-ray emission as a function of source power through the use of multiwavelength 
observations. 



I. INTRODUCTION 

A central question that is been debated in the Fermi 
era regards the location of the blazar 7-ray emission 
site: is the 7 - ray emission of powerful blazars pro- 
duced inside the sub-pc size broad line region (BLR) 
or further out at scales of ~ 1 - few pc (fig. 1) where 
the IR photon field of the dusty molecular torus (MT) 
dominates over that of the UV field of the BLR? In 
the first case the 7 - ray emission is most probably ex- 
ternal Compton (EC) scattering of the ~ 10 eV BLR 
photons [IH] , while in the second the seed photons for 
the EC emission are the ~ 0.1 eV photons emitted by 
the dust [9] in the MT [8]. The issue of the energy dis- 
sipation location is connected to the jet formation and 
coUimation process and, as we propose, can be 

approached through a SED diagnostic we introduce. 



A. The near camp: Few hour GeV variability 
puts the blazar inside the BLR. 

Fermi has detected flares with decay time down to 
~ 3 hours, comparable to the telescope sky scan- 
ning period. Such variations have been seen in 3C 
454.3 H [33], PKS 1454-354 [J, 3C 273 [3], PKS 
1502+106 pi, PKS B1222-H216 [II]. The short vari- 
ability timescale has been used to argue that the emis- 
sion is produced within the BLR (e.g. [33]): For a jet 
with r = 5 = 10, where F is the bulk Lorentz factor of 
the flow and S the corresponding Doppler factor, vari- 
ability times of tyar = 10^ s (~ 3 hours), correspond 
to a maximum source size of r = ctyar^ = 3 x 10^^ 
cm. Assuming a jet half-angle Oj^-t = 0.17/F [T7], an 
upper limit on the distance of the blazar from the 
central engine is i? = r/9jet = rF/O.l? = 3.5 x 10^"^ 
cm. This distance is comparable to the BLR size 
Rblr ~ 1 — 3x 10^^ cm[Tni and, therefore, it is plausi- 
ble that the GeV emission is produced inside the BLR. 
Essentially this argument is based on the assumption 
that the entire cross section of the jet is emitting. If 
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FIG. 1: Is the blazar GeV emission produced inside the 
sub-pc broad line region or further out at distances ~ 1 
- few pc, comparable to the size of the molecular torus? 
This is an ongoing debate with strong cases presented by 
both camps. Here we introduce a diagnostic based on the 
spectral energy distribution (SED) and we show how it 
can be used to locate the blazar emission site. 



however one allows smaller parts of the jet to produce 
the 7-ray emission (e.g. [14]) such short variations can 
take place further out. 



B. The far camp: Optical/VLBI polarization 
and GeV/X-ray data put the blazar few pc from 
the central engine. 

The sub-pc scale energy dissipation is challenged 
by observations that put the emission at few pc dis- 
tance from the black hole (e.g. [H [5l [H] [24l ) . In 
several cases, optical polarimetry during an optical - 
7-ray flare showed a polarization behavior similar to 
that observed in simultaneous high frequency VLBI 
coming from the 43 GHz core, several pc away from 
the central engine. Because of the similar optical and 
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VLBI polarization behavior, the optical emission is 
constrained to emerge from the VLBI core and be- 
cause the optical and 7-ray variations are seen to be 
simultaneous, the 7-ray emission is also constrained 
to emerge from the VLBI core at a distance of few to 
several pc from the black hole. Additional arguments 
placing some of the 7-ray flares at ^ 10 pc have been 
advanced for 3C 454.3 [521 [SO]- These arguments are 
based on the similar behavior of optical and millimeter 
light curves and explain the GeV emission as EC scat- 
tering of photons emitted from the dust of the MT. 
Although at distances of ^ 10 pc from the central 
engine variability events are expected to be of the or- 
der of ~ 10 days, significantly shorter variations (day 
long or even down to a few hours) are seen (e.g. PKS 
1510-089 [2S]). These can be explained if the observed 
variations come from a fraction of the jet cross section. 

Recently, a diagnostic has been proposed [TT] for the 
location of the blazar GeV emission that is based on 
the energy dependence of Fermi flux variations. This 
diagnostic makes use of the fact that if the blazar emis- 
sion takes place inside the BRL the cooling is done on 
the ~ 10 eV BLR photons, while if it takes place out- 
side the BLR, the most abundant seed photons are 
IR photons (~ 0.1 eV of the MT. In the first case 
cooling takes place at the onset of the Klein-Nishina 
cross section and GeV variability is achromatic, while 
in the second case cooling takes place in the Thom- 
son regime and variability is faster at higher energies. 
Although this is a powerful diagnostic, it can only be 
applied to a small number of bright Fermi flares. 

II. AN SED-BASED 7-RAY LOCATION 
DIAGNOSTIC AND ITS APPLICATION 

In powerful blazars, the emission consists of two 
spectral components. The low frequency one peaks 
at Vg ~ 10^^ Hz [inj and is attributed to synchrotron 
radiation; the high frequency one peaks at ~ 10^^ 
Hz J£J and is usually attributed to inverse Compton 
emission from external photons that are produced in 
the broad line region or in the MT. 

If F is the bulk Lorentz factor and S the usual 
Doppler factor of the jet flow, and 7b is the Lorentz 
factor of the electrons responsible for the synchrotron 
and EC SED peaks, then, 

= #-7?<5, (1) 



ec = ^eo7fc'5^ (2) 

where es is the synchrotron peak energy, ec is the EC 
peak energy, eo is the characteristic energy of the ex- 
ternal photon field, all in units of the electron rest 
mass, B is the magnetic field permeating the emission 



region and Bcr = (mlc^) / (eh) — 4.4 x 10"^^ G is the 
critical magnetic field. Note that equation ^ is valid 
only if the scattering of electrons with Lorentz factor 
7b takes place in the Thomson regime. The condi- 
tion for this, eo7b(5 ^ 1 can be written with the help 
of equation ^ as (ecEo)^^^ ^ 1- The highest possi- 
ble energy external seed photons are UV line photons 
with eo ~ 10~^, which means that the scattering is in- 
deed in the Thomson regime as long as ec ^ e,^"^ « 10^. 
This corresponds to an energy of « 5 GeV or ~ 10^''' 
Hz. Given that in most cases powerful blazars peak 
at lower i/c, with < i^c >~ 10^^ Hz [Tn|, the scattering 
of the electrons producing the EC peak is well within 
the Thomson regime. Taking the ratio of the two peak 
energies we obtain: 

B _ Ae^esBcr /„v 

The same ratio B /S can be obtained from the expres- 
sion for the Compton dominance fc, the ratio Lc/Ls 
of EC to synchrotron luminosity: 

where Uq is the external photon field energy density in 
the galaxy frame, Uq = (4/3)J7or^ is the external pho- 
ton field energy density in the jet comoving frame and 
Ub ~ B^/{8tt) is the magnetic field energy density 
[121 US]- Solving equation Q for and equating to 
equation ([s]), we obtain our final expression 

U^ = = 3.2 X 10^ G (5) 

eo V 67r ec i^c,22 

where fci is the Compton dominance in units of 10, 
Vc,22 is i^c in units of 10^^ Hz, and 1^3,13 is Vs in units 
of 10^'^ Hz. JVote that the RHS contains only observ- 
ables. It is the RHS that informs us about the ratio 
of the square root of the energy density over the peak 
energy of the seed photons (seed factor, SF, in Gauss) 
available at the location of the j-ray emission.lt is this 
information that we can use for understanding where 
the emission comes from. 



A. The SF in the BLR and in the MT. 

Reverberation mapping finds that Rblr. ~ 1 — 
3 X lO^^Ly^g cm where L^; 45 is the accretion disk 
luminosity in units of 10^^ erg s~^, and that a 
fraction ^ ~ 0.1 of is reprocessed by the BLR 
[3 OS]. The energy density for the BLR is then 
Uo = ^Ld/(47ri?2c) w 0.29 - 2.6 x 10"^ erg cm'^ 

^ 1 /2 

[15] . Because R oc Lj , Uq is the same for sources 
of different luminosities. The BLR SED in the galaxy 
frame can be approximated by a blackbody peaking 
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FIG. 2: I/c as a function of the SF for a number of sources, 
color-coded for redshift, for which Vs, Vc and k axe found 
from their SEDs. The solid line is the average from |16| . 
the areas marked by dotted lines correspond to the 7-ray 
emission site being inside the BLR (left) or the MT (right). 



at I'o ~ l-Si^LyQ (^o ~ 3 x 10 ^, [32 ). Using these we 
obtain SF - 1.8 - 5.5 x 10^ G. 

Because of the larger distance of the MT from the 
central engine, reverberation mapping has only been 
performed for Seyferts (e.g. [31]), lower luminosity 
and therefore smaller MT size sources. These stud- 
ies, along with IR interferometric studies (e.g. j20| ) 
are also in agreement with an i? cx L^J"^ scaling, 
suggesting that Uq is the same for sources of differ- 
ent luminosities. Adopting the results of [23 for the 
blazar 4C 21.35, a blackbody of temperature T = 1200 
K (eo = 5.7 X 10-'') and L = 7.9 x lO'*^ erg s'^ 

1 — 2 pc, we obtain 
With these we ob- 



emittcd from a radius of 
Ua = 0.5 - 2.2 X lO'^ergs cm 
tain SF ~ 4.1 — 8.2 x lO'' G. A comparison of these 
indicative ranges is shown in Figures 2 and 3. 



B. Can SSC be dominant? 

In the previous discussion we considered only ex- 
ternal photons. Here we briefly discuss under what 
conditions synchrotron photons would be the domi- 
nant seed photons. A lower limit for he synchrotron 
photon density Ug in the comoving frame (assuming 

(5 = r is 



47rc3C,r6' 



(6) 



where t^ar is the observed variability timescale and Lg 
the luminosity of the synchrotron component. For this 



to dominate over UqT^ , the comoving external photon 
field energy density, we require 



r < 



(- 



12.6 



L 



sA7 



(7) 



For our adopted BLR range, C/q « 0.29— 2.6 x 10 ^ erg 
cm~^, this is equivalent to F < 8.4 — 11.0 ^iVra/e/t 
Similarly for our adopted MT range, Uq ~ 0.5 — 
2.2 X 10~'^erg cm^'' this is equivalent to F < 11.4 — 
13.7 -'^s ^Vra/e/i- VLBI studies of superluminal speeds 
in FSRQs [e.g. see figure 24 of [17] show that for most 
FSRQs, F > 10, with values reaching up to 40. This 
raises the possibility that our diagnostic is relevant for 
sources that have relatively high F, and that for the 
slower sources SSC may dominate. 



III. 



PRELIMINARY RESULTS. 



Some preliminary results using the sample of |26| 
with detections in the 2FGL or BAT are shown in Fig- 
ures 2 and 3. The average values of k '-^few, Vg — IQ^^ 
Hz, and Vc = 10^^ found for powerful blazars [TO] 
correspond to a SF (solid line) that is closer to the 
MT range than that of the BLR. The fact that the 
two bands corresponding to the SF we estimated for 
the BLR and the MT occupy the central part of the 
observed SF range is encouraging and indicates that 
FSRQs experience a seed photon environment that is 
within the range of what is expected from our simple 
BLR and MT considerations. We note that estimat- 
ing Vg and Vc requires good multi-wavelength cover- 
age, and errors of a factor up to ~ 10 are possible, 
particularly in the current estimates of v^- This sug- 
gests that the actual range of SF is significantly more 
narrow than the current data suggest, likely giving 
closer agreement with the range denoted by the bands 
in Figures 2 and 3 once we have completed our SED 
analysis on very well sampled sources. 

In future work, we plan to use all available data, 
including archival data from sources like NED with 
newer catalogues from WISE, Planck, SWIFT, BAT, 
LAT to produce SEDs with simultaneous or contem- 
poraneous data for as many sources as possible and 
derive their Vg , and k following |26j , along with er- 
ror estimates for these quantities. We anticipate that 
this work will significantly expand the sample in |26j 
which was restricted to coverage of the synchrotron 
component only. 

To estimate the actual range of the SF in the BLR 

and the MT, we need to know how the quantity 
1/2 

C/q /cq is expected to vary radially. This requires 
information on the BLR and the MT stratification, 
something that is not yet well constrained [Tl fT^E^ ) 
yet. For the modeling component of the project, we 
plan to explore how the SF range widens if one con- 
siders the stratification of different line photons in the 



eConf ClllllO 



4 



2011 Fermi & Jansky: Our Evolving Understanding of AGN, St Michaels, MD, Nov. 10-12 



51 - 









' blr' 


'mt' 






I *I 








✓ 

/ A 






A ^ 

✓ 








1 * 1 




AA 


1 1 A 








■ 






1 1 A 








• 




1 1 










• ; • ? • 












• i 








»l 1 










V- !. i 












z = 1 .5-2 


^ • 1 1 




• z>2 


✓ 




A BATFSRQ 









BLR (e.g. P7, and the fact that powerful IR compo- 
nents with T down to ^ 300 K have been detected in 
radio-loud AGN 122j. 
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FIG. 3: Same as Figure 2, but for sources with z > 2. 
There is a clear correlation between the total IC power and 
the seed factor SF, moving from values consistent with the 
BLR to the MT at highest energies. 
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